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CAN NEUTRAL AND IONIZED POLYCYCLIC AROMATIC HYDROCARBONS BE CARRIERS OF THE 
ULTRAVIOLET EXTINCTION BUMP AND THE DIFFUSE INTERSTELLAR BANDS? 
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ABSTRACT 

Up to now, no laboratory-based study has investigated polycyclic aromatic hydrocarbon (PAH) 
species as potential carriers of both the diffuse interstellar bands (DIBs) and the 2175 A UV bump. 
We examined the proposed correlation between these two features by applying experimental and 
theoretical techniques on two specific medium-sized/large PAHs (dibenzorubicene C30H14 and hexa- 
benzocoronene C42H18) in their neutral and cationic states. It was already shown that mixtures of 
sufficiently large, neutral PAHs can partly or even completely account for the UV bump. We in- 
vestigated how the absorption bands are altered upon ionization of these molecules by interstellar 
UV photons. The experimental studies presented here were realized by performing matrix isolation 
spectroscopy with subsequent far-UV irradiation. The main effects were found to be a broadening of 
the absorption bands in the UV combined with slight redshifts. The position of the complete tt - tt* 
absorption structure around 217.5 nm, however, remains more or less unchanged which could explain 
the observed position invariancc of the interstellar bump for different lines of sight. This favors the 
assignment of this feature to the interstellar PAH population. As far as the DIBs arc concerned, 
neither our investigations nor the laboratory studies carried out by other research groups support a 
possible connection with this class of molecules. Instead, there are reasonable arguments that neutral 
and singly ionized cationic PAHs cannot be made responsible for the DIBs. 

Subject headings: astrochemistry — dust, extinction — ISM: molecules — methods: laboratory — 
molecular data — techniques: spectroscopic 



1. INTRODUCTION 

Mid-infrared (MIR) emission bands related to aromatic 
C-H and C-C vibrational modes were found in several 
astrophysical environments, ranging from planetary neb- 
ula, reflection nebulosities, circumstellar disks to even ac- 
tive g alactic nuclei (see, e.g., the recent review of lTielen^ 
|2008[ ). Inferred from their presence, polycyclic aromatic 
hydrocarbons (PAHs) are believed to be an important 
component of interstellar dust, and, after H2 and CO, 
they arc probably the most abund ant molecules in the in - 
terstellar medium (ISM; see, e.g., iLeger fc PugetlllQSl . 
Although these MIR bands do not permit an identifica- 
tion of specific PAH molecules, general information, such 
as size and charge state distributi on or the aliphatic - 
to-aromatic ratio, can be inferred (jDraine fc Lil l2007h . 
Based on energetic arguments, the MIR emission bands 
were attributed to aromatic mo lecules contain ing about 
(or at least) 50 carbon atoms (jTielensI I2008D . Besides 
the elaborately investigated IR properties of PAHs, it 
is equally important to perform experimental studies in 
other wavelength regimes to test whether the PAH model 
is compatible with observations. Usually, neutral or ion- 
ized PAHs possess no or only weak permanent dipole 
moments. Therefore, strong spectroscopic fingerprints 
due to pure rotational transitions in the radio range arc 
not expected. On the other hand, electronic tt - tt* tran- 
sitions give rise to distinct absorption bands in the UV- 
VIS. The ISM displays two clear spectroscopic features in 



this wavelength range: a growing number of yet unidenti- 
fied diffuse interstellar bands (DIBs) above 400 nm and a 
broad bump at 217.5 nm with almost fixed position, but 
varying strength and width for different lines of sight. 
These features are superimposed on a background with 
increas i ng ex tinction at shorter wavelengths (see, e.g., 
lDrainel[200l . Below 400 nm, the extinction curve is 
smooth and no narrow bands that can be r elated to poly- 
atomic ga s phase m olecules we re found (iClavton et ahl 
[2003; Gre del et al.i r2011: SalamaeTaD 120111). Unfortu- 
nately, electronic spectra of astrophysically relevant large 
PAHs and mixtures of them measured in the labora- 
tory under appropriate conditions are almost unavail- 
able, mainly because suffici ent quantities of su i table sam- 
ples are difficult t o obta in. iRuiterkamp et al.] ()2002D and 
iHalasinski et al.l ()2003[) studied the UV-VIS absorption 
properties of a few selected matrix-isolated PAHs con- 
sisting of max. 48 C atoms. These studies were guided 
by the search for possible DIB carriers and were biased 
toward more exotic or elongated structures because of 
their strong absorption bands in the VIS. By measuring 
the absorbancc of mixtures o f hot gaseous PAHs with 
mean sizes around 31 C atoms. [Joblin et aTl (|1992D found 
that the strongest spectral feature is a broad bump at 210 
nm, very close and similar in shape to the interstellar ex- 
tinction bump at 217.5 nm. Recently, we obtained anal- 
ogous results for mixtures of comparably dimensioned 
PAHs which were isolated in solid Ne at low tempera- 
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ture (jSteglich et al.ll2010D . A strong UV bump was found 
around 200 nm accompanied by a rather featureless ex- 
tinction curve if molecules with less than 22 C atoms were 
absent from the mixture. Considering semi-empirical cal- 
culations, we argued that the bump position may shift to 
the red, eventually reaching 217.5 nm, for species large 
enough to be associated with the MIR emission bands 
{> 50 C atoms). This idea is furthermore supported by 
calculations applyi ng time-dependent de nsity functional 
theory (TD-DFT). iMalloci et all ()2004[ ) showed that a 
PAH mixture, comprising a limited number of different 
structures, could also give rise to a bump around 217.5 
nm. 

While more than 400 DIBs are known today, none 
of their carriers which are thought to be large carbon- 
bearing molecules in the gas phase, such as PAHs and 
their ions , have been identified yet (for a review on 
DIBs see ISarri I2006D . An unambiguous identification 
requires the comparison with experimental spectra of 
the species under consideration. Usually, this demands 
for the application of elaborate experimental techniques, 
like cavity ring down spectroscopy carried out at low 
temperature in the gas phase. So far, the spectra of 
only a few small PAHs consisting of less than 25 C 
atoms h ave been recor ded with this method (for a re- 
view see iSalamal l2008f ) . When in their neutral charge 
state, such small molecules display their strongest ab- 
sorption bands in the UV. As already mentioned, no 
sharp absorption bands were found in the UV part of 
the interstellar extinction curve i mplying a very small 
abundance of sma l l neutral PAHs (iClavton et al.l 12^0031 : 
IGredel et al.llMlt ISalama et al.l[201lD . The cations, on 
the other hand, show additional strong bands in the VIS 
wavelength range. However, these bands are usually too 
broad compared to the widths of most DIBs. Recent at- 
tempts to search for bands of the naphthalene (CioHg) 
and anthracene (C14H10) cations i n interstellar spec- 
tra were proven to be u nsuccessful (jGalazutdinov et al.l 
[20Tlt ISearles et al.llMl . 

This article presents new experimental results on the 
UV-VIS absorption properties of medium-sized to large 
PAHs. Supported by theoretical calculations, we at- 
tempt to solidify the connection between large PAHs and 
the UV bump. Bearing in mind that, for different lines 
of sight, the molecular population is affected by differ- 
ent physical conditions, like the strength of the UV field, 
it is important to address the question whether ioniza- 
tion has an effect on the bump position of PAH mix- 
tures. We approach this issue experimentally through 
photoionization of two matrix-isolated PAHs, dibenzo- 
rubiccne C30H14 (DBR) and hexabenzocoronene C42H18 
(HBC), which display strong absorption bands around 
217.5 nm in their neutral state. For the first time, com- 
plete electronic spectra of PAH ions spanning the full 
UV-VIS range down to 200 nm will be presented. The 
spectrum of the HBC cation will be of p articular interest 
as, in addition to its superior stability ()Trov fc Schmidt! 
I2006D . this molecule has a size comparable to the PAH 
size that was inferred from the aromatic IR emission 
bands. 

2. PREDICTIONS INFERRED FROM TIME-DEPENDENT 
DENSITY FUNCTIONAL THEORY 

2.1. UV Bump at 217.5 nm 



In this section, we shortly address the problem of the 
UV signature of PAHs (the bump below 230 nm) and 
its size dependence using computational methods. Pre- 
viously, we already calculated the absorbance of arti- 
ficial PAH mixtures with differe nt mean sizes by ap - 
plying a semi-empirical method (jSteglich et al.l l2010( ) . 
The low computational costs of this method, however, 
may be accompanied by rather unprecise band positions 
and strengths. Furthermore, size-dependent effects could 
partially be related to the model itself. This can be 
avoided by using more elaborate ab initio and TD-DFT 
methods which, on the other hand, can only be ap- 
plied to a limited number of molecules. IMalloci et al.l 
()2004D calculated the absorption spectra of 40 PAHs 
in different charge states using a special treatment of 
TD-DFT implemented in the Octopus soft ware package 
(jCastro et al.ll2006HMarques fc Grossl2004l ). These spec- 
tr a. were later collected in an o nline database and used 
bv lCecchi-Pestellini et al.l (|2008f ) to simulate artificial ex- 
tinction curves for a dust model that includes PAHs. 
In this model, the UV bump and the non-linear far-UV 
(FUV) rise of the interstellar extinction curve were solely 
attributed to electronic transitions in PAHs. 

Here, we used the TD-DFT implern entation of the 
Gaussian09 software (iFrisch et al.ll2009l ) which relies on 
basis sets to create the molecular orbitals. The computa- 
tional effort of this method scales steeply with the num- 
ber of electronic transitions to be calculated and, there- 
fore, it is much more expensive compared to the Octopus 
implementation when the spectrum has to be obtained 
over a large energy range. At present, it is the most accu- 
rate method to calculate electronic spectra of large PAHs 
down to wavelengths of ca. 200 nm, i.e., for energies 
smaller than 6.2 eV (5 iJ.m~^). We focused on seven se- 
lected PAHs with compact structures, possibly represent- 
ing the average PAH population of the ISM. As indicated 
by astronomical observations in the mid-IR, interstellar 
PAHs are rather compact than e longated (iLeger et al.l 

19891) The B3LYP functional (iStephens et al.l 119941: 

Beckel[T993l ) in conjunction with the 6-31-|-G(d) basis set 
was used to optimize the ground state geometries and 
to calculate the vertical electronic transition energies af- 
terward. Additional vibrational excitations cannot be 
treated in this approach. However, at higher energies, 
a substantial band broadening is expected, resulting in 
spectra without discernible vibrational pattern even for 
cold matrix-isolated or gas phase molecules (see, e.g.. 
Section [3]) . 

Figure [1] displays the calculated electronic absorption 
spectra of compact PAHs belonging to the Dgh (top 
panel) and D2h (bottom panel) point groups. It should 
be taken into account that the uncertainty for the band 
positions at this level of theory can typ ically reach up 
to 0.3 eV (see, e.g., iHirata et al.l Il999( ). At the mo- 
ment, we limit the discussion to neutral species, as the 
ions demand higher computational efforts due to their 
open-shell structure. For each series (D2h and Dgh), the 
same symmetry-related selection rules apply, allowing di- 
rect comparison among the differently sized species. The 
theoretical stick spectra (oscillator strength vs. energy, 
not shown) have been convoluted with Lorentzian func- 
tions with full width half maximum (FWHM) of 3000 
cm~^ to obtain the spectra shown by solid lines. Similar 
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absorption spectra can be expected for weakly bonded 
PAH clusters where vibrational patterns are smoothed 
out due to the van der Waals interactions among the 
molecules. To illustrate this point, the absorption spec- 
tra of thin films of coronene and HBC are shown for 
comparison (dash-dotted lines). For the spectra repre- 
sented by dashed lines, the FWHM for each transition 
was chosen to be proportional to the number of energet- 
ically lower-lying transitions while keeping the area of 
each band proportional to the corresponding oscillator 
strength. The purpose of this convolution is to reflect 
the increased lifetime broadening at higher energies and 
to visually enhance the weaker bands at longer wave- 
lengths. This also resembles more closely the situation 
of matrix-isolated molecules as will be seen in Section 
[3l However, one should consider that the band intensity 
(integrated cross section) of real molecules in the visi- 
ble is usually distributed over several vibrational bands 
reducing the peak intensities. 

Looking at both PAH series, a strong UV peak around 
217.5 nm for molecules with more than sa 32 C atoms 
(ovalenc) is ap parent which is in agreement with the re- 
sults of Mallo ci et al.l (|2004l ) who showed that the cations 
and anions also display this UV peak (see also Fig. [2]). 
From the computational point of view, the transitions 
of large PAHs that produce the UV bands around 217.5 
nm are almost not affected when one electron is removed 
from the highest occupied orbital or added to the low- 
est unoccupied orbital. Differences in the UV spectra of 
different charge states should mainly concern the band- 
widths. This problem cannot be assessed with theo- 
retical calculations and will be addressed in Section [31 
While the feature around 217.5 nm seems to be a com- 
mon property among different sufficiently large PAHs, 
bands at longer wavelengths depend more strongly on 
the specific molecular structure. The averaging over sev- 
eral PAH molecules will produce a rather smooth extinc- 
tion curve with a broad UV bump, possibly around 217.5 
nm, as main feature. Even though the species we chose 
for the calculations are not large enough (due to compu- 
tational limitations), a shift of their UV bump to even 
longer wavelengths, i.e., beyond 220 nm, can be expected 
for increasing molecular sizes. It is known that also 
amorphous graphitic material and nanoparticles, usu- 
ally containing both aromatic and aliphatic components, 
can feature a similar UV bump in the wavelength range 
between 220 and 260 nm (Sc hnaite r ct al. 1996, 1998!; 
iLlamas-Jansa et al.l l2007| ) . The varying band position 
might, at least partly, be related to the average size of 
the aromatic planes. However, the width of the bump 
in these nanoscopic and macroscopic materials is often 
quite broad due to several possible effects, such as parti- 
cle agglomeration, van der Waals and covalcnt bonding 
between the aromatic planes, or just broad size distribu- 
tions of the aromatic planes. 

2.2. FUV Rise 

In the wavelength range below 200 nm, the spectra of 
large PAHs share further characteristics. These charac- 
teristics can only be obtained by theory, as almost no ex- 
perimental data are available. For the calculation of the 
FUV absorption bands of PAHs, the method described 
in the previous section is computationally too expensive 
due to the high density of states in this energy regime. 
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Fig. 1. — Theoretical electronic absorption spectra of different 
compact PAHs with Dsh (top panel) and D2h (bottom panel) sym- 
metry calculated at the B3LYP/6-31-|-G(d) level of theory. A con- 
volution using Lorcntzians with energy-dependent widths has been 
applied to visually enhance the weaker bands at longer wavelengths 
(dashed lines). For comparison, the absorption spectra of thin films 
of coronene and HBC arc shown (dash-dotted lines). The position 
of the interstellar UV bump is marked at 217.5 nm, plus-minus a 
reasonable 0.3 eV uncertainty accounting for the inaccuracy of the 
TD-DFT calculation. 



The determination of electronic transitions should now 
be based on the T D-DFT implementa t ion of the Octopus 
software package () Castro et al.l 120061 : iMaroues fc Gross 

m Compared to the calculations bv iMalloci et~ 
, we chose slightly different parameters. The Oc- 
topus code relies on numerical meshes instead of basis 
sets to calculate the wave functions in real-space and 
real-time. Followed by an initial very short electromag- 
netic pulse, exciting all frequencies of the system, the 
time-dependent dipole moment is calculated from which 
the linear optical absorption spectrum, comprising all 
possible electronic bound-bound transitions, can be de- 
rived. The crucial parameters of the calculation are the 
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functional (here B3LYP), the volume of the numerical 
box (sphere around each atom with radius 3 A), the grid 
spacing (0.25 A), the time integration length (20 ?i/eV), 
and the time step (0.002 h/eV). 

Here, we only want to emphasize the most important 
and securest information that can be extracted from 
these calculations and that concerns the compatibility 
of the FUV absorption properties of large PAHs with 
the interstellar extinction. A mo re extensive discussion 
can be found in the publication bv lCecchi-Pestellini et al.l 
(I2OOI . In this paper, different FUV rises that are ob- 
served for the interstellar extinction in different lines of 
sight were interpreted to stem from different ionization 
states in the same PAH population. The calculated PAH 
spectra showed a tendency for a reduction of the absorp- 
tion cross section in the gap between the highest-energy 
TT - TT* transitions and the a - a* FUV rise when going 
from anions over neutral molecules to cations. However, 
this effect seems to be more pronounced in smaller PAHs 
and the majority of the 40 different molecules in this 
model contained less than 32 C atomsQ Also, no physi- 
cal explanation was given for this behavior. 

Using the Octopus code, we calculated the spectra of 
neutral, cationic, and anionic HBC which are displayed 
in Figure pi The spectr a are v ery similar to the ones ob- 
tained bv lMalloci et al.l (|2004[ ) who used the LDA (local- 
density approximation) functional. The widths of the 
absorption bands are purely artificial and depend solely 
on the integration length used in the calculation. The 
area of each band, however, is directly related to the os- 
cillator strength of the corresponding transition(s). The 
shown spectra of HBC below 8 //m~^ are dominated by 
TT - TT* transitions. The calculated energies of their band 
positions, however, are somewhat too low as becomes ob- 
vious upon comparison with the measured spectrum of a 
film-like HBC deposit. Also, the calculated strength of 
the transition at 2.3 /im~^ is underestimated by a factor 
of two. Note that the broad bands of the HBC film will 
become much narrower if the molecules are isolated from 
each other in the gas phase or in a low-temperature ma- 
trix. By comparing with experimental gas phase spectra 
of anthracene, it was shown that the general trend of the 
broad a - a* hump above 8 /xm~^ as calculated with 
the Oc topus code is in go od agreement with measure- 
ments (jMalloci et al.ll2004[) . In the energy range above 2 
/im~^, the calculated spectra of differently charged HBC 
molecules are principally very similar to each other as 
this range is dominated by transitions where electrons 
are promoted from occupied to unoccupied molecular 
orbitals, i.e., electrons overleap the band gap (B- and 
C-type transitions). Further to the red (< 2 /im^^), ad- 
ditional transitions to the semi-occupied orbital of the 
cation (A- type transitions), as well as combined transi- 
tions to and from the semi-occupied orbital of the anion, 
give rise to additional bands which are absent for the 
neutral molecule. 

For the interstellar extinction, no correlation could be 
observed between the bump at 217.5 nm and the non- 
linea r FUV rise at shorter wavelengths (jCardelli et al.l 
Il989( ). In order to ascribe the UV bump to a PAH fea- 

^ For details, see the online database at 
littp: / /astrochemistry.ca. astro. it/database/pahs. html. 
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Fig. 2.— TD-DFT (B3LYP) spectra of neutral and ionized 
hexabenzocoronene (HBC) calculated wit h the Octopus code 
HCastro et al.|[2006l : IMargues fc Gross|[200l ) for the elucidation of 
the FUV properties of large PAHs. The shaded area roughly marks 
the absorption gap between the tt - tt* bands and the onset of 
the (7 - fT* hump. Note that the energetic positions of the tt - 
TT* resonances below 6 \j,ra~^ are actually underestimated by this 
method which becomes apparent in com parison with the me asured 
spectrum of a film-like deposit of HBC llSteglich ct al.|[2010l '). The 
mean interstellar extinction curves for two different reddening pa- 
rameters i?v (wh ere Ry = 3.1 is the in terstellar average) are shown 
for comparison l lCardelli et al.|[l98g|) . The jy-axes are labeled in 
units of cross section per H atom cth or C atom ctc, respectively. 

ture, the FUV absorption properties of PAHs have to 
be compatible with the observed FUV rise of interstel- 
lar dust which is weakest for lines of sight with strong 
reddening (see Figure [2]) . As is apparent from the cal- 
culated HBC spectra, the onset of the a - a* hump is 
at rather high energies (8 ^m^^). The absorption gap 
between the strongest tt - tt* resonance and the a ~ g* 
absorption onset fits quite well to the observed FUV ex- 
tinction properties of the interstellar dust. Furthermore, 
the theoretical spectra illust rate that the effect which was 
described bv [Cecchi-Pcstell ini et al.l ([2008), i.e., the dif- 
ferently pronounced gap for different charge states, does 
not necessarily have to be true for larger species. Cer- 
tainly, alternative scenarios that could explain the ob- 
served variations of the strength of the non-linear FUV 
rise, such as changes in the composition of other dust 
components, especially aliphatic components in carbona- 
ceous grains, have to be explored. 

3. LABORATORY-MEASURED SPECTRA OF SELECTED 
NEUTRAL AND IONIZED PAHS 

3.1. Matrix Isolation Spectroscopy 

This section addresses the question how photoioniza- 
tion influences the UV absorption bands of large PAHs, 
especially the bump around 217.5 nm. According to the 
ca lculations presented before and the theoretical studies 
bv lMalloci et al.l (|2004[ ) the energies and strengths of the 
short-wavelength B- and C-type transitions of large neu- 
tral PAHs (A < 450 nm in the case of HBC) are almost 
unaffected upon removal or attachment of one electron. 
Potential differences, e.g., concerning the band shapes 
and positions, however, are only accessible via labora- 
tory measurements. For this purpose, we used matrix 
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isolation spectroscopy (MIS) to measure the wavelength- 
dependent photoabsorption between 190 and 850 nm of 
HBC and DBR molecules isolated in solid Ne at 7 K. The 
number of Ne atoms relative to the number of molecules 
to be investigated (isolation ratio) was on the order of 
15000 for DBR and 30000 for HBC, respectively. After 
preparation of the matrix doped with the neutral PAH 
molecules, the sample was exposed to the FUV radia- 
tion (10.2 - 11. 8 eV) of a hydrogen- flow discharge lamp 
(|Warneckl[T96l to create ions from the neutral parent 
molecules. During the irradiation of typically 15 min, 
the FUV intensity on the sample surface was roughly 
IQis _ iQie photons m'^ s'^ The applied low FUV 
doses, as well as the inability of the molecules to dif- 
fuse through the matrix exclude the possibility that more 
complex chemical alterations than simple one-photon in- 
duced reactions (here only ionization) take place. The 
setup for MIS a nd FUV processing i s described in more 
detail elsewhere (jSteglich et al.ll201ll ). Parts of the spec- 
trum of ionized HBC in Ar matrix, which have been cho- 
sen for presentation because of a better si gnal-to-noise ra- 
tio, w ere recorded with a different setup ( Bouwman et al.l 
|2009( ). Due to the close proximity of the molecules in 
the matrix, recombination reactions between already ion- 
ized molecules and released electrons limit the maximum 
ion yield to values somewhere below 20% after a cer- 
tain FUV dose is reached. To disentangle the absorp- 
tion bands of ionic and neutral species which, especially 
in the UV range, overlap each other, the following pro- 
cedure was performed (see Fig. |3l). The transmission 
spectrum measured after irradiation was first baseline- 
corrected against the spectrum before irradiation, i.e., 
the spectrum of the neutral molecules. The spectrum 
thus derived in units of absorbance (red curve in Fig. [3]) 
displays positive features from species created during the 
irradiation (ions) and negative features from species that 
were destroyed (neutrals) . Then, the contributions of the 
destroyed neutral molecules were simply subtracted by 
using an accurate fit obtained from the complete spec- 
trum before irradiation. After subtraction, the spectrum 
contains only bands from species that were created dur- 
ing the irradiation. Actually, this procedure works quite 
well in the UV range despite the strong overlap of bands 
of the neutral and ionized molecules. For instance, even 
the sharp lines at 308 nm (3.25 /im~^) and 283 nm (3.53 
/xm~^) from the OH radical (green curve in Fig. [3]) which 
was created by photodissociation of trace amounts of wa- 
ter are properly reproduced. Problems can only occur in 
spectral ranges where the spectrum of the neutral pre- 
cursor contains sharp bands that are represented by only 
one or two data points. For DBR, this is the case only 
at 540 nm (1.85 fim~^). 

The electrons released from the neutral parent 
molecules are usually trapped at defects and impu- 
rities in the matrix. Therefore, cations are gener- 
ally created in exc ess compared to anions (see, e.g., 
iSalama fc Allamand ola 1991, 1991). The ratio between 
cation and anion number densities depends on several 
parameters, such as the purity and crystallinity of the 
matrix, the isolation ratio, or the electron affinities of 
the molecules and impurities. The possible contribution 
of anions to our spectra will be discussed in the following 
subsections. 
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Fig. 3. — Procedure for obtaining the full-range spectrum of ion- 
ized PAH molecules, here DBR, in an inert gas matrix is illustrated 
(see the text for an explanation). The spectrum of neutral DBR 
shown at the top was baseline-corrected to account for light scat- 
tering from the matrix. The finally obtained spectrum of ionized 
DBR at the bottom displays an uncorrected interference pattern 
in the VIS. This has been removed in the spectrum shown in Fig. 

m 



3.2. Dibenzorubicene (C^qHh) 

The measured spectra of neutral and ionized matrix- 
isolated DBR are displayed in the top panel of Figure 
m Detailed information on the spectroscopy, both in ab- 
sorptio n and emission, of n eutral DBR is given in another 
paper (jRouille et al.|[20Tll ). DBR can be regarded as a 
planar fragment of the fuUerene C70 (with H-saturatcd 
dangling bonds), which was recently dis covered in the 
planetary nebula Tc 1 (jCami et al.H2010[ ). The neutral 
molecule has several strong absorption peaks in the UV, 
one specific at 212 nm in close proximity of the 217.5 nm 
interstellar bump. Therefore, it is an ideal test molecule 
to investigate how the high-energy tt - tt* transitions are 
affected by ionization. 

The theoretical spectra of neutral, cationic, and an- 
ionic DBR are presented in the bottom panel of Figure 
[31 They were obtained in the same manner as the dashed 
spectra of Figure [H i.e., using the TD-DFT approach of 
the Gaussian09 program, applying the 6-311-|--|-G(2d,p) 
basis set, and performing a convolution with Lorcntzians, 
where the widths were proportional to the number of 
energetically lower-lying states. The optimized ground 
state geometry of DBR in all three charge states was 
found to be C2h- The electronic ground state symmetries 
are "'^Ag (neutral), ^Bg (cation), and ^Au (anion). For 
the neutral and the cation, electronic transitions from 
the ground state to excited states of Au and Bu sym- 
metry are dipole-allowed. For wavelengths shorter than 
500 nm, the electronic excitation schemes of the various 
transitions are very similar for the neutral and the cation. 
Computational differences mainly arise from the broader 
cation bands which are caused by the convolution proce- 
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dure and the fact that the cation has several low-energy 
transitions unlike the neutral molecule. Above 500 nm 
the theoretical spectrum of the DBR cation is dominated 
by weak transitions where electrons arc promoted to the 
scmi-occupicd orbital, e.g., D5(Bu)^Do(Bg) at 673 nm, 
D6(Bu)^Do(Bg) at 613 nm, D8(Bu)^Do(Bg) at 567 nm, 
and Dio(Bu)<— Do(Bg) at 540 nm. Two further dipole- 
allowed (Bu) and two dipole-forbidden (Ag) transitions 
are predicted beyond 900 nm, outside of the accessible 
scanning range. The symmetry-related selection rules 
that apply to the anion are somewhat different. Due to 
its ungerade ground state, dipole-allowed transitions are 
only possible to Ag and Bg states. Therefore, its calcu- 
lated spectrum deviates more strongly from the spectra 
of the other two species. In view of the experimental 
spectrum of ionized DBR and the lack of strong absorp- 
tion bands beyond 600 nm, we conclude that much less 
anions were formed than cations and we ascribe the mea- 
sured spectrum solely to the DBR cation. 

Except for a sometimes wrong ordering of states, it 
can be seen that the applied TD-DFT method works 
fairly well in the case of the neutral molecule. The 
Si(Bu)<— So(Ag) transition, for instance, is predicted at 
536 nm, while the origin band is actually measured at 
540 nm. Two close-lying bands, S2(Bu)<— So(Ag) and 
S4(Bii)^So(Ag), are calculated to be at 448 and 427 nm. 
The weaker one at 448 nm (S2), with observed origin at 
434 nm, is responsible for the vibrational pattern on top 
of the broad band (S4) that can be seen in the experimen- 
tal spectrum between 360 and 440 nm. These two bands 
have their counterparts with similar band positions and 
strengths in the cationic molecule, where they appear 
broader and slightly red-shifted in the measurement. To 
the red of these two bands, two weak and broad features 
(marked with * and **) are observed around 620 and 690 
nm which could correspond to the Dio-s— Dg and Dg^— Do 
transitions. 

In the wavelength range between 340 and 200 nm, more 
than 100 single electronic transitions arc predicted by 
TD-DFT for the cation. Roughly 45 of them have non- 
zero oscillator strengths and actually contribute to the 
observed broad tt - tt* structure peaking at ca. 215 nm. 
Except for minor peak shifts, this band system resembles 
a broadened version of the analogous tt - tt* structure of 
the neutral precursor which essentially can be explained 
by similar excitation schemes of electronic C-typc tran- 
sitions for both molecules. Because of a stronger broad- 
ening for all UV bands, the bump peaking at 215 nm 
appears even more pronounced for the cation compared 
to the neutral. 

3.3. Hexabenzocoronene (C42H1S) 

The spectroscopy of matrix -isolated neutral H BC is 
discussed in detail elsewhere (jRouille et al.l 120091) . Be- 
cause of its all-benzenoid structure, it is believed that 
HBC is more abun dant in space compar ed to other 
similarly-sized PAHs ()Trov fc Schmidtll2006[) . The calcu- 
lated spectrum of neutral HBC, obtained at the B3LYP 
/ 6-31-|-G(d) level of theory, reveals a good match 
with the MIS measurements (see Figure [S]), both in 
terms of relative band strengths and positions. Neutral 
HBC has a Dgh ground state geometry. The intensity 
of its S4(Eiu)^So(Aig) transition, calculated to be at 
364 nm, is actually distributed over a complicated vi- 
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Fig. 4. — Top: measured spectra of matrix-isolated neutral 
and ionized dibenzorubicene (C30H14, DBR) at low temperature. 
Sharp features from the OH radical created by the dissociation of 
trace amounts of water are marked accordingly. The position of 
the interstellar UV bump is indicated by the vertical dashed line. 
Bottom: calculated spectra of DBR {B3LYP/6-311-f-(-G(2d,p)) in 
different charge states. 

brational pattern that involves aig modes of vibration 
and vibrational excitations belonging to an energetically 
lower Biu state, to which electronic transitions from 
the ground state are dipolc-forbiddcn. Also the weak 
Si(B2u)'<— So(Aig) transition, for which a peak of an un- 
defined vibrationally excited state has been measured at 
434 nm, borrows some intensity from the S4(Eiu) state. 
(Actually, the Eiu s tate could also be S3. For details, see 
iRouille et al.l ()2009[ ).) Quite well reproduced by the cal- 
culation, further strong bands arc located around 217.5 
nm. 

The spectroscopic interpretation of the spectrum of 
ionized HBC in solid Nc is more complicated. The degen- 
erate ground states of the HBC cation and anion cause 
Jahn- Teller interaction, which effectively leads to a re- 
duction of the ground state geometry from Dgh to D2h. 
This also complicates the assignment of measured ab- 
sorption bands on the basis of DFT calculations as these 
methods can already fail to correctly predict the dis- 
torted geometry of a degenerate ground sta te (see, e.g., 
iDavidson fc BordenI [19831 iRuss et aL|[200l . The com- 
parison of the measured spectrum of ionized HBC with 
the TD-DFT calculated spectra (bottom panel in Fig- 
ure [S]) does not permit to exclude the presence of anion 
bands, in contrast to the case for DBR. The band sys- 
tem with observed origin around 830 nm in the Ne matrix 
could have its theoretical counterpart in either cation or 
anion bands predicted around 740 nm. Using the real- 
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time, real-space implementation of TD-DFT (see Figure 
[2]), however, these anion bands are calculated to be much 
further to the red around 1100 nm (0.88 /im"^), while 
the cation bands remain around 760 nm (1.31 ^m~^). As 
already stated, cation bands usually dominate the spec- 
tra of matrix-isolated photoionized PAHs. Therefore, we 
tentatively attribute all observed bands of the photopro- 
cessed matrix (except the OH lines) to the HBC cation 
and discuss the results accordingly. 

The electronic ground state of the HBC cation trans- 
forms according to the Big irreducible representation of 
the D2h point groupQ Electronic transitions to ^Au, 
^B2u, and ^B^u states are dipole-allowed. The theoret- 
ical spectrum contains several A-type transitions to the 
semi-occupied orbital at wavelengths longer than 500 nm. 
The lowest Di(B2g)^Do(Big) transition is actually cal- 
culated to be far beyond 2 fim. The observed band pat- 
tern between 750 and 835 nm may originate from the fol- 
lowing calculated transitions: D6(B3u)^Do(Big) at 732 
nm with oscillator strength / = 0.038, D7(Au)<— Do(Big) 
at 724 nm with / = 0.039, and D8(B3u)^Do(Big) 
at 686 nm with / — 0.081. Another transition 
Dg(Aii)-(— Do(Big) is predicted at 542 nm, close to the 
measured broad band with maximum at 531 nm. How- 
ever, its strength (/ = 0.025), if it can be ascribed to the 
measured band, is considerably underestimated which, 
e.g., could be related to an underestimation of the Jahn- 
Teller distortion in the ground state by DPT. 

At wavelengths shorter than 400 nm, the density of 
electronic transitions rises quickly. A detailed spectro- 
scopic analysis is hardly practicable anymore. The com- 
parison of ionized HBC in two different matrix environ- 
ments (Ar and Ne) illustrates an important aspect con- 
cerning the absorption spectra of PAH cations in this 
wavelength range. As it was the case for DBR, the spec- 
trum of ionized HBC contains broad bands which resem- 
ble broadened versions of analogue bands found in neu- 
tral HBC. Especially, the tt - tt* structure around 217.5 
nm is even more apparent for the ionized molecule. If the 
strong band broadening could be ascribed to a matrix- 
induced effect, larger bandwidths for the Ar compared to 
the Ne environment would be expected, as interaction- 
related effects are usually more pronounced for atoms 
with higher polarizabilities. The shapes and widths of 
these features, however, do not seem to depend on the 
rare gas atoms surrounding the molecule (for A < 550 
nm), indicating an intrinsic effect, e.g., very short life- 
times of the excited electronic states. On that sup- 
position, low-temperature gas phase spectra of ionized 
medium-sized to large PAHs, such as DBR and HBC, at 
wavelengths shorter than approximately 400 nm would 
feature band shapes and widths similar to what has been 
measured for molecules isolated in the inert gas matrix. 
However, a small wavelength shift, usually to the red, 
has to be considered in that case. The matrix-induced 
redshift becomes obvious upon comparison of the Ne and 
Ar measurements. We will take a closer look at this in 
the following subsection. 



^ All given symmetries refer to the molecule being oriented in 
the y-z-plane (x=0). The z-axis traverses the C atoms and defines 
the direction in which the molecular structure is elongated. 
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Fig. 5. — Top: measured spectra of matrix-isolated neutral and 
ionized hexabcnzocoroncnc (C42Hig, HBC) at low temperature. 
Sharp features from the OH radical created by the dissociation 
of trace amounts of water are marked accordingly. The position 
of the interstellar UV bump is indicated by the vertical dashed 
line. Bottom: calculated spectra of HBC (B3LYP/6-3H-G(d)) in 
different charge states. 



3.4. Implications for the Carriers of the Diffuse 
Interstellar Bands 

A connection between large PAHs and the DIBs is of- 
ten assumed. Nevertheless, DBR can be excluded as 
possible DIB carrier. While the case o f the neutral 
mole cule is discussed in another article (jRouille et al.l 
|2011[ ). cationic DBR does not exhibit sufficiently nar- 
row bands - at least for wavelengths shorter than 900 
nm. The first electronic transition of HBC in the gas 
phase causes some narrow and very weak bands between 
410 and 435 n m (24390 - 22990 cm ~^) not able to repro- 
duce any DIB (jKokkin et al.|[2008l ). Stronger absorption 
bands appear further in the UV where no DIBs were ob- 
served. The HBC cation, on the other hand, displays 
a system of rather narrow bands between 750 and 850 
nm. In Ne matrix, the strongest features were found 
at 830.8 nm (12037 cm"!) and 827.3 nm (12088 cm"!). 
The availability of the spectra of ionized HBC in matrices 
with different polarizabilities allows us to perform a lin- 
ear extrapolation procedure to obtain the band positions 
in the gas phase. This pro cedure is described in another 
article ( Gredel et al.ll20l"lT ). The wavenumber position of 
a certain band in the gas phase I'gas can be approximated 
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11800 12000 12200 12400 12600 12800 13000 13200 13400 



^gas 



VAr 



(1) 



where VMe,Ar are the respective band positions in the 
matrices. The polarizabiHties amount to aAr = 1-64 A'^ 
and aNe^ 0.397 A^ (jRadzig fc SmirnovllMl . The ex- 
trapolated gas phase positions and approximate oscilla- 
tor strengths for bands of ionized HBC are summarized in 
Table [H For the bands between 750 and 850 nm (13400 - 
11800 cm~^), the positions have been obtained by fitting 
the experimental spectra to Lorentzian functions. For 
the broad asymmetric band at 531 nm (18840 cm~^), 
only the maximum position has been taken. We did not 
measure the absorption cross sections directly in the MIS 
experiments because the results are very often not accu- 
rate. Instead, the approximate oscillator strengths were 
determined by comparing the integrated band intensity 
with the integrated intensity of the strong band system 
between 200 and 260 nm. The accumulated oscillator 
strength of this band system was taken from the TD- 
DFT calculation (/ ?a 7). The uncertainties of the deter- 
mined /-values mainly reflect the expected inaccuracies 
of the calculated transition strengths. The inaccuracies 
were obtained for the neutral molecule upon comparison 
with experimental data of HBC in a solvent. 

The linear extrapolation procedure on the basis of the 
absorption spectra is furthermore illustrated in Figs. [5] 
and [T] The bands that are displayed in Fig. [6] may 
actually originate from more than one electronic tran- 
sition which could also explain small differences in the 
matrix-induced wavelength shifts (see the third column 
in Table [T|). Moreover, these bands seem to be superim- 
posed on a very broad band with a tail extending down 
to 600 nm (16670 cm~^; see also Fig. [5]). The spectra 
of Fig. [7] probably display a vibronic progression with 
rather broad bands. The strongest feature is an asym- 
metric band peaking at 531 nm (18840 cm~^) in Ne. 
Its band shape and bandwidth (k, 300 cm"-'^) are iden- 
tical in both matrix environments and are expected to 
be the same in the gas phase. A comparable band shape 
has been found a t 462 nm in Ar matrix for the smaller 
coronene cation (jSzczepanski fc Valal I1993D , which has 
an electronic structure similar to the HBC cation. 

The comparison of the band positions that were ex- 
trapolat ed for the gas phase w i th the synthetic DIB spec- 
trum of iJenniskens fc DesertI (|1994f FI demonstrates that 
ionized HBC can safely be excluded as a potential DIB 
carrier. Basically, the obtained gas phase positions and 
/-values could be used to search for interstellar features 
of HBC cations, but we are skeptical whether such an 
endeavor would be successful. A search for interstellar 
HB C in its neutral cha rge state was already unsuccess- 
ful (jCredel et al.ll201l|) . We want to stress again that 
because of its size and superior stability (most stable 
PAH with 42 C atoms) one would expect to observe fea- 
tures of neutral or charged HBC on the interstellar ex- 
tinction curve if the PAH hypothesis, i.e., the assignment 
of the DIBs to large PAHs, would be correct. However, 
in view of the expected chemical variety for PAHs of this 

^ The DIB spectrum has been calculated with the 
most recent data from the online catalog available at 
http://leomd.arc.nasa.gov/DIBcatalog.html. 
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Fig. 6. — Comparison of absorption spectra of ionized HBC in Ar 
and Ne matrices with the sp ectrum of the diffuse interstellar bands 
llJenniskens fc Desert) [T994I ). Presumably, all measured bands are 
due to the HBC cation. The gas phase positions of the strongest 
HBC features which were obtained by a linear extrapolation proce- 
dure are indicated (see also Table [TJ. The fitted Lorcntzians were 
only used to determine the band positions. 

size and the failure to detect HBC in interstellar spec- 
tra it seems unlikely that a few selected PAH molecules 
(out of several million possible structures) could be abun- 
dant enough to be responsible for one or more of the 
roughly 400 DIBs. Additionally, PAH spectra can be 
very complex with complicated vibrational and interme- 
diate level structures which do not seem to be present 
in interstellar spectra. This is especially true for PAHs 
with molecular structures that are less symmetric than 
what is usually studied. Furthermore, electronic bands of 
cations appearing in the VIS can be very broa d (see, e.g., 
iBiennier et al][200l lSukhorukov"eFaI|[200l making it 
impossible to connect them with the narrower DIBs. In 
summary, we would expect that the population of large 
interstellar PAHs produces a smooth extinction curve as 
it was already observed for mixtures o f sufficiently large 
PAHs in inert gas matrices (jSteglich et al..,201(]h . 

4. CONCLUSIONS 

We presented experimental and computational spec- 
troscopic studies of electronic transitions in the UV-VIS 
of PAHs comprising up to 42 C atoms (66 C atoms in the 
calculations). Principally, these species share two com- 
mon features in the UV. The first one is a strong tt - 
TT* resonance below 230 nm, which will be particularly 
pronounced for a mixture of (probably rather compact) 
PAHs and whose exact position will depend on the molec- 
ular mean size. It will probably peak around 217.5 nm if 
the size of the PAHs is such that they can be connected 
to the aromatic MIR emission band s. As was demon- 
strat e d for natural PAH mixtures (jEhrenfreund et al.l 
[1991 IStegTich et al.l[2010l) . narrow bands below 400 nm 
appear on the collective tt - tt* slope of the absorp- 
tion curve if sufficient amounts of small molecules (< 
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TABLE 1 

Extrapolated Gas Phase Positions and Approximate Oscillator Strengths for Bands of Ionized HBC (see also Figs. [61&0 



Band Position in Matrix [cm i'Ng — i^Ar Extrapolated Gas Phase Approx. Osc. 



Ar (12 K) 


Ne (7 K) 


[cm-i] 


Position [cm ^] 


Strength (+20%) 


11889 ± 9 


12037 ± 8 
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12084 ± 14 


0.0065 


11935 ± 9 


12088 ± 7 


153 
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0.0032 
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Fig. 7. — Comparison of absorption spectra of the HBC cation 
observed at 536 and 531 nm in Ar and Ne matrices with the spec- 
trum of the diffuse interstellar bands ( Jenniskens & Desert 1994). 
The gas phase position (529 nm) was determined by a linear ex- 
trapolation. 



22 C atoms) are present. However, sensitive searches 
failed to detect such sharp UV absorption features in 
the diffuse IS M (IClavton et al.l [20031 : IGredel et al.l[20TTI : 
[ski ania ct al, 20 lit ) indicating the absence or low abun- 
dance, respectively, of free-flying small PAHs which is 
in accordance with the size constraints of interstellar 
PAHs as inferred from the MIR observations. Strin- 
gent constraints on the abundance of spe cific PAHs have 
also been obtained by other studies (,Kokkin et al.1120081 : 



iRouille et al1[200l [Pilieri et al.l[2009l l. At shorter wave- 
lengths, the spectra are dominated by the broad a - 
a* hump whose properties are quite common, especially 
among larger molecules. Its absorption onset is around 
125 nm (8 /.tm"^) as was demonstrated for HBC. The 
gap between the highest-energy tt - tt* feature and the 
a - <j* onset fits quite well to the observed interstellar 
extinction curve, even for lines of sight displaying a low 
FUV extinction. 

Using HBC and DBR as examples, we investigated ex- 
perimentally how photoionization influences the tt - tt* 
absorption bands for wavelengths shorter than 400 nm. 
It was found that the main effect is a broadening [3 in 
combination with slight wavelength shifts of the individ- 
ual bands. The position of the complete tt - tt* struc- 
ture around 217.5 nm remains rather unaffected which 
can be explained with similar excitation schemes of C- 
type transitions in this energy range for the neutrals and 
the cations. The spectra of ionized HBC in different in- 
ert gas matrices illustrate another important point. The 
widths and shapes of the broad bands at wavelengths 
shorter than 600 nm (in the case of HBC) do not de- 
pend on the specific matrix environment. These bands 
are much broader than what is expected from typical 
matrix-induced broadening effects indicating very short 
lifetimes of the excited states which are not entirely 
caused by the interaction with the matrix atoms. (Al- 
ternatively, other broadening mechanisms, like Franck- 
Con don vibronic broade ning, may be considered (see, 
e.g.. ISassara et alll2001[) .) Therefore, absorption spec- 
tra of cold ionized HBC in the gas phase (below 600 
nm) will most likely resemble the presented MIS spectra. 
However, a small shift of the band positions to shorter 
wavelengths can be expected. 

* The issue of the band broadening of the first electronic transi- 
tions and non-radiative relaxation rate s in small neutral and ion- 
ized PAHs was recently summarized by IPino et al.l 1 120111) . 
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If the interstellar UV bump at 217.5 nm is indeed a 
collective feature of the interstellar PAHs that are ob- 
served through their IR emission signatures, its observa- 
tional constraints could be explained as follows. The 
varying width for different lines of sight might be an 
indication for different degrees of ionization. Alterna- 
tively, this could just hint at different size distributions 
of the molecules. (Also, both explanations could be pos- 
sible.) Its seemingly fixed position can probably be ex- 
plained by a more or less stable mean size of the PAHs 
in the mixture. Within the observational and computa- 
tional limitations, it was already shown that the inter- 
stellar abundance of PAH s could suffice to explain the 
observed buinp str ength (jCecchi-Pestellini et al.l 120081 : 
iSteglich et al.ll20Tol) . As suggested by astronomical ob- 
servations, the outflows of carbon-rich stars injec t hydro- 
genat ed amorphous carbon (HAC) into the ISM (jMathisI 
119901) . Recent laboratory investigations of such ma- 
terials indicate the formation of more aromatic com- 
pounds under long-term FUV light exposure and the 
appearanc e of a bump-like abs orption feature around 
217.5 nm (IGadallah et al.ir2011f ). However, this feature 
is yet too broad in comparison with the interstellar UV 
bump, probably because the electronic tt resonances in 
the HACs are still affected by chemical and physical 
bonding between the aromatic subunits. Furthermore, 
these laboratory materials seem to need slightly more 
carbon than available for interstellar dust in order to re- 
produce the observed strength of the interstellar feature 
(jGadallah et al.l l20 lit ) which may be rationalized by a 
high content of aliphatic structures. Further chemical 
alteration by, e.g., FUV irradiation or interstellar shocks 
in the ISM could lead to the destruction of such struc- 



tures and the release of aromatic units, i.e., large PAHs, 
with sizes sufficient to explain the MIR emission as well 
as the UV bump and the smooth extinction curve below 
400 nm. These large PAHs, either neutral or ionized, 
with variable hydrogen saturation, may then represent a 
more or less stable end product in the ISM. 

Concerning the PAH hypothesis, we found no evidence 
that neutral and singly ionized cationic PAHs studied so 
far by us or other laboratory groups could be responsible 
for any of the DIBs. Especially, the case of HBC, 
which is rather large and by comparison very stable, is 
compelling. The apparent absence of electronic features 
of neutral and cationic HBC in the diffuse ISM suggests 
that, due to molecular diversity, the abundances of 
individual PAH species of this size are probably too 
low to allow detection in the UV-VIS spectral domain. 
Therefore, we conclude that the DIBs (at least the 
prominent ones) are most likely not signatures of the 
interstellar PAH family, if we consider only neutral and 
singly ionized cationic species with normal hydrogen 
saturation. 
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